We studied the turbulence spectrum of the local interstel- 
lar plasma in the direction of PSR J0437-4715, on the basis of 
our observations and those reported earlier by others. We com- 
bine these data to form a structure function for the variations of 
phase along the line of sight to the pulsar For observations that 
did not report them, we infer modulation indices from a theoret- 
ical model. We find that all of the observations fit a power-law 
spectrum of turbulence with index n - 3.46 ± 0.20. We suggest 
that differences among reported values for scintillation band- 
width and timescale for this pulsar arise from differences in 
observing parameters. We suggest that refractive effects domi- 
nate for this line of sight, with refraction angle about twice the 
. . . , diffraction angle at 330 MHz observing frequency. We suggest 
that the scattering of this pulsar lies in a layer of enhanced tur- 
bulence, about 10 pc from the Sun. We propose that the flux 
' variations of the extragalactic source PKS 0405-385 arise in 
. ' flie same scattering layer 
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Abstract. 



1. Introduction 

Study of the interstellar medium (ISM) for nearby pulsars leads 
to understanding of the local interstellar medium. As Shishov 
et al. (I2003> showed, analysis of multi-frequency observations 
of pulsar scintillation is critical to understanding the turbulence 
spectrum of the interstellar plasma. Multifrequency observa- 
tions are able to detect details that cannot be seen from ob- 
servations at a single frequency. For example, construction of 
the phase structure function from multi-frequency observations 
of PSR B0329H-54 resulted in the detection of strong angular 
refraction along the line of sight to this pulsar. This structure 
function also showed that the power-law index of the spectrum, 
3.50 + 0.05, differs significantly from the Kolmogorov value of 
11/3. The Kolmogorov spectrum describes much data on pul- 
sar scintillation quite well, over a large range of spatial scales 
(Armstrong et al. 1995 Shishov & Smirnova 2002 1, but in par- 
ticular directions, and particularly along short lines of sight 
which sample only a small part of the interstellar medium, the 
spectrum can differ from the Kolmogorov form. 

In this paper we study the turbulence spectrum of plasma 
along the line of sight to PSR J0437-4715. This is one of the 
closest pulsars. It has distance R = 150 pc, and its transverse 
velocity is 100 km/sec (van Straten et al. .2001. 1. It is quite 
strong over a wide range of observing frequencies. We com- 
bine the data in the accompanying paper (Gwinn et al. 2005 
hereafter Paper 1; also Hirano 2001 1 with observations by oth- 
ers, and construct the phase structure function for this pulsar 
in the time and frequency domains. From examination of the 
structure functions, we conclude that refractive effects are im- 
portant for this line of sight, with refraction angles about twice 
the diffraction angle. Based on comparison of data for the flux 
variation of the extragalactic source PKS 0405-385 with the 
scintillation parameters of PSR J0437-4715, we suggest that 
that the intensity variations for both are caused by the same 
scattering layer, located near the Earth at a distance ^ 10 pc. 



2. Observational Data 

2.1. Scintillation Parameters 

Observers commonly describe a pulsar's scintillation by its 
characteristic scales: scintillation bandwidth and scintillation 
timescale. These describe the decline of the autocorrelation 
function of the intensity, with frequency and time lag, from 
its peak at zero time and frequency lag. Observers usually 
normalize the autocorrelation function by its peak value. This 
is the square of the modulation index, if the data are nor- 
malized by the mean intensity (/). The modulation index is: 
m = V^/^'^^'W^/(^)- Scintillation timescale tiss is the time 
lag where the autocorrelation falls to 1/e of this central value, 
and scintillation bandwidth Aviss is the frequency lag where 
it falls to 1 /2 of that value. For fully sampled data, in strong 
scattering where differences among wave paths are many ra- 
dians, the modulation index is m - I. If the data do not span 
several scintillation bandwidths, or frequency scales, then the 
modulation index is less than 1, as discussed in the Appendix 
below. 

Although the autocorrelation contains a great deal of infor- 
mation on the structure function of the density inhomogeneities 
responsible for scattering (Shishov et al. 2003 1, the characteris- 
tic scales of scintillation carry a very limited part of that infor- 
mation. A more detailed study of the structure function of the 
inhomogeneities requires the modulation index and the form of 
the autocoiTelation function. 

2.2. Summary of Observations for PSR J0437-4715 

Tabled compares measurement of scintillation timescale and 
bandwidth for PSR J0435-4715, from the Uterature. For easier 
comparison, we scaled these measurements to a single observ- 
ing frequency / - 330 MHz, using the scaling relations ap- 
propriate for a Kolmogorov spectrum, n - 11/3: tjss /''^ 
and Avjss Fig- Q shows the same data, in graph- 

ical form. Column 7 of Table [l] gives the observing band- 
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width scaled to 330 MHz, corresponding to the horizontal lines 
shown in Fig. [J The measurements fall into two clear groups: 
those showing narrow-band scintillations with Av ~ I MHz, 
and those showing broader-band scintillation with Av ^ 15 
MHz. Observations cannot determine scintillation bandwidth if 
greater than the observing bandwidth B. Fig.^shows this up- 
per limit for each observation. As the figure shows, Issur ( 2000 1 
and Paper 1 had sufficient bandwidth to detect the wider-band 
scintillations; other observations could detect only the narrow- 
band scintillations. 

The scintillation bandwidth of the narrow-band scintillation 
varies between epochs, as is found for some other nearby pul- 
sars (Bhat et al. 1999 1. For PSR J0437-4715 this variabiHty 
is most pronounced in the observations of Gothoskar & Gupta 
J2000t . who found that the measured Aviss spans a range of 
more than an order of magnitude within 2 days. Paper 1 finds a 
value for the narrow-band scintillation at both epochs that lies 
within this range; Gothoskar et al. observed between their two 
epochs. Nicastro & Johnston J1995> and Johnston et al. ( 1998J 
observed variation of only 30% over more than 2 years. The 
ranges of values the two groups measure for the narrower Avjss 
overlap, and have similar centroids. We will show in ^3.3l and 
l^that all observed scales of scintillation can be explained by 
one power-law spectrum of inhomogeneities. 



In the difftactive model, decorrelation with frequency arises 
from changes in the scattering angle with frequency. The fre- 
quency difference Afj corresponding to the scale b is 

A/rf ^ c{kbflR oc b^f, (4) 

where l/{kb) is the typical scattering angle 9sc, k is the wave 
number and R is an eff'ective distance to the scattering layer 
(Ostashev & Shishov [D771 Shishov et al. .2003J . Thus, for 
diffractive scintillation, one obtains: 

mh) = (fo/ffmf)- (5) 

In a model including strong angular refraction, another re- 
lation between Af and b is realized. We introduce the typical 
refractive angle Ore/- If dref ^ ^likb), the displacement of a 
beam path with changing frequency converts the spatial pattern 
to structure in frequency. The frequency difference, Afr, corre- 
sponding to the scale b is then (Shishov et al. .2003j : 

Afr^c{kb)m,f)^bf. (6) 

So for the case of strong angular refraction: 

A/.(/o) = (/o//)'A/(/). (7) 

We will use these relations in our analysis of the data. 



3. Structure Function 

3.1. Construction of the Structure Function 

Only a few of the observations of PSR J0437 -47 15 provide the 
shape of the correlation function, which can be used to con- 
struct the phase structure function Ds{t, f). As Shishov et al. 
J200 3^ showed, Ds{At) can be obtained for small time lags Af 
from the correlation function B/(Af) of intensity variations: 

B/(0) - BAAt) 

D,(Af) = ^^-^^yyj^ forAf<f„5, (1) 

where tiss is the characteristic time scale of variations caused 
by electron density inhomogeneities in the ISM. In the fre- 
quency domain we have: 

B/(0) - BAAf) 

Ds(Af)^ '^jyj^ ^' forAf<Av,ss, (2) 

where Aviss is the characteristic frequency scale of scintilla- 
tion, and Af is a frequency difference. 

To scale all data to a single reference frequency /o, we must 
rescale the structure functions at frequency / appropriately: 

DAfo, Ar(/o), A/(/o)) = At, Af){fo/ff. (3) 

The time difference, Af, characterizes the same spatial scale in- 
dependently of observing frequency, At(fi)) = Af(/). However, 
^.fif) A/(/())- This is because wave scattering by inhomo- 
geneities with a spatial scale b determines the decorrelation of 
intensity fluctuations at frequency difference Af. Moreover, the 
relation between b and Af depends on observing frequency /. 
One may consider two different models for the relation between 
b and Af: diffractive and refractive (Shishov et al. i2003i) . 



3.2. Inferred Scintillation Parameters of PSR 
J0437-4715 

As shown in detail in Table ^ scintillation parameters for 
PSR J0437-4715 obtained at different frequencies and differ- 
ent times differ very strongly even when referred to a single 
frequency: by about a factor of 30. We observe that there can 
be at least two explanations for such behavior One is that mul- 
tiple scales of scintillation are present; one might expect such 
a model in the interstellar Levy flight proposed by Boldyrev & 
Gwinn J2003t . for example. In this case we might see fine fre- 
quency structure caused by occasional large-angle deflections, 
leading to large propagation times, as well as the coarser struc- 
ture originating from more typical deflections. 

An alternative explanation is that we do not observe the 
fundamental scales of scintillation: for example, if the band- 
width of the receiver is less than the actual diffractive scale, 
then one observes only the tail of the fast variations of intensity 
in frequency and time domains. This tail has modulation index 
m « I. Usually observers do not record the modulation index, 
but rather only the scintillation bandwidth and timescale as de- 
scribed above. When the bandwidth of the analysis is less than 
the diffractive scale of scintillation, Av/ss, these scales do not 
correspond to the actual time or frequency scale of scintillation. 

3.3. Inferred Modulation Index 

Only two observations of those Usted in Table^had sufficiently 
wide observed bandwidth to detect the fundamental scale of 
scintillation: Issur (2000) and Paper 1. All other observers had 
bandwidths several times less than the characteristic frequency 
scale of scintillation. Unfortunately they did not publish ob- 
served values of the modulation index, but we can determine 
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the expected modulation index theoretically if we know the ra- 
tio of the observing bandwidth B to the scintillation bandwidth 
Avjs s , and the ratio of the time span of a scan T to the scintil- 
lation timescale f/j 5 . 

As is shown in the Appendix, the modulation index is ex- 
pected to be defined by the largest of the ratios of these param- 
eters, if both ratios are less than one. We note that the modula- 
tion index will be different for diffractive and refractive models 
(see Appendix). In the presence of strong angular refraction, 
when the refraction angle is much more than the scattering an- 
gle, 9ref » Odif, a refractive model is required (Shishov et al. 

We determined the expected modulation indices from the 
observing bandwidth and time span, and show the results in 
Table 12 For all of the observations, T/tjss > 1, but for some 
we have B/Avjss < 1. In this case, the expected modulation 
index is given by (see Appendix): 

(B/Av/ss)*""^^^^ diffractive model, (8) 



{n(n + 2)) 
0.7 



refractive model. 



(9) 



(n(n - 1))J 

Here, n is the power-law index of the spectrum of density fluc- 
tuations in the ISM: for example, for the Kolmogorov theory 
n = 11/3. In the standard calculation of 2D (frequency and 
time) autocorrelation functions of scintillation spectra, and de- 
termination of tiss and Avjss from them, the mean intensity 
from the observations is subtracted from the intensity data be- 
fore autocorrelation (Cordes[l986j. Equation |8] or |9] holds in 
this case. Table |2l includes the expected value of the modula- 
tion index for both models, m^rf and ;«£,,-. These values were 
calculated from the expected values of the scintillation band- 
width Av° and the timescale t^^^ . We found these expected val- 
ues from the values reported by Paper 1: Avjss = 16 MHz and 
tjss - 17 minutes for / - 328 MHz. We scaled these values 
to the observing frequency / using the scaling laws appropri- 
ate for a Kolmogorov spectrum, n = 11/3: Av° oc and 
tiss °^ /'^- For the refractive model, the scaling law is actually 
Av" oc Z"*^, but this does not change our estimates by more than 
10%. 

We note that the frequency scale measured at 152 MHz 
agrees well with that extrapolated from 328 MHz, but the es- 
timate of the decoiTelation time has a large eiTor bar because 
the ratio T/tjss is small. In this case, it is better to use tjss 
defined at the level of 0.5, rather than 1/e. We did not include 
the values of tjss and Avjss measured at 660 MHz by Johnston 
et al. 



I1998> because the ratio B/Avjss is very small for those 
measurements, so that the mean intensity cannot be defined re- 
liably. 

Paper 1 reports two frequency scales: Avi = 16 MHz and 
Av2 - 0.5 MHz, both at an observing frequency of f - 328 
MHz. The small-scale scintillation was detectable with a re- 
duced observing bandwidth of B = 8 MHz, as we discuss fur- 
ther below. They found that the modulation index for this struc- 
ture is m 0.2, which is close to what we would estimate from 
Eq.|9l Gothoskar & Gupta (20001 report Av^^ * 2.5 MHz at 
the same frequency. This differs markedly from the value re- 
ported by Paper 1 . Their receiver bandwidth was less than half 



the expected decorrelation bandwidth, but they had finer fre- 
quency resolution than Paper 1 . 

The structure function of intensity fluctuations with fre- 
quency coiTesponds to the actual structure function of interstel- 
lar scintillation only when the frequency bandwidth of the ob- 
servations, B, is significantly larger than the frequency scale of 
the scintillations. In that case one obtains the true characteristic 
frequency scale. In the case of limited observing bandwidth B, 
one observes only a tail of the actual frequency structure func- 
tion, corresponding to the small scales of variations. The scales 
of these variations obtained from the autocoiTelation function 
will depends on the parameters of the observation, and can be 
used for construction of a composite structure function only by 
using the proper modulation index (see Ea. llOt . In the case of 
insufficient bandwidth and time span of observation, the scin- 
tillation bandwidth and scintillation timescale obtained from 
the autocorrelation function can differ strongly from the actual 
ones. 

4. Results 

We constructed a composite time and frequency structure func- 
tion using coiTelation functions at different frequencies, and 
converting them to a frequency /o, using Eqs. Ill3l and either 
|8lor|9l and either|5lor0for the frequency correlation function. 
We used the frequency correlation function obtained by Issur at 
152 MHz, frequency and time correlation functions at / = 327 
MHz from the paper of Gothoskar & Gupta ( 20-00) ^ and data at 
328 MHz of Paper 1 . We normalized the correlation functions 
by dividing them by their values at zero lag, so all are normal- 
ized at zero lag. We did so because we do not know the actual 
values of mean and rms intensity for the different observations. 
In this case. 



D,(Af,A/) = (1 - B{At,Af))m^ 



(10) 



We chose to use /o - 1000 MHz, to simplify comparison with 
structure functions for other pulsars determined by Shishov et 
al. (EOOD and Smirnova et al. (2006). 

In Fig. 12] we present the phase structure function. Note 
that the data from Gothoskar & Gupta (2000.1 and Paper 1 
agree very well. The expected modulation index, as calculated 
in Table 121 has been used in setting the overall offset of the 
data of Gothoskar & Gupta. We have estimated the value of 
the structure function at 436 MHz based on the characteristic 
timescale reported by Johnston et al. (1998), and associating 
that point with the expected modulation index given in Table 
121 for the refractive model. This point is in good agreement 
with the other points, and that good agreement leads us to the 
conclusion that all spatial scales of inhomogeneities probed in 
these observations arise from a single spatial spectrum. A fit to 
all points gives us the slope a - \ .46+0.20, which corresponds 
to n = 3 .46 + 0.20. To convert this structure function of time lag 
At to a function of spatial scale (, we use the simple conversion 
as ^ = Af Viss, where Viss is the speed of the diffraction pat- 
tern relative to the observer. To choose Vjss correctly, we must 
know the location of the scattering material along the line of 
sight. We consider this location and the appropriate conversion 
m Obelow. 
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In Fig. 13 we show the frequency structure functions for two 
models: a model with a strong angular refraction (upper panel) 
and a diffractive model (lower panel). The data have been com- 
bined using Eqs.^-I^] All of the data agree well for both mod- 
els, but the slope of the structure function (a = 1.44 + 0.03) 
is about the same as for the temporal structure function. This 
leads to the conclusions that we have a refractive model: for 
the diffractive model the slope should rather be half that of the 
temporal structure function (Shishov et al. l2003> . 

5. Discussion 

Several components of the ISM are responsible for the inter- 
stellar scintillation of pulsars and extragalactic radio sources. 
Component A is homogeneously distributed up to a distance of 
order 1 kpc from the Sun., and is distributed in the space outside 
the spiral arms of the Galaxy (Cordes 1985 Cordes et al. 1991, 
Pynzar' & Shishov I1997> . At larger distances component B, 
located in the spiral arms of the Galaxy, gives the primary con- 
tribution to scintillation effects. A third component lies about 
10 pc from the Sun, where it contributes an enhanced level of 
turbulence (Jauncey et al. 120001 Dennett-Thorpe & de Bruyn 
120021 Rickett et al. 2002). We refer to this component as com- 
ponent C. Component C is primarily responsible for the scintil- 
lation of extragalactic radio sources, because the large angular 
sizes of the sources suppress the influence of other components. 
For the closest pulsars, the scintillation effects caused by com- 
ponents A and C can be comparable, and we must investigate 
especially the relative contribution of these two components in 
scintillation, for each pulsar. 

Pulsar PSR J0437-4715 is one of the closest pulsars. 
Estimation of the relative contributions of components A and 
C for this object depends on the relative values of the spatial 
diffractive scale b and the Fresnel scale bfr - (R/^y^^, where 
k - 2n/Ais the wavenumber, and R is the distance to the scat- 
tering layer for the pulsar. 

5.1. Model for Scattering by Component A 

Suppose first that component A provides the primary contribu- 
tion to the scintillation of PSR J0437-4715. Using the distance 
to the pulsar of L = 150 pc (van Straten 2001 1, we obtain the 
value of the Fresnel scale for component A (here L ^ R) at 
frequency / = 330 MHz, and find bp, = (R/k)^^^ ^ 10" cm. 
Using the pulsar velocity of V = 100 km/s (van Straten UOOU . 
we obtain the characteristic spatial scale of the diffraction pat- 
tern bdif - Vtiss ~ 10'" cm. Because bf,- > bdif, the scintilla- 
tion must be saturated. Using the value of bdif we can estimate 
the scattering angle O^c ~ 1/ikbdif) ~ 0.3 mas. The expression 
for the scintillation bandwidth can be obtained from Equation 

il 

Av%f ^ c/(R0i). (11) 

Using Eq[n]and the scattering angle Osc - 0.3 mas, we obtain 
the estimated scintillation bandwidth for the diffraction pattern: 
Av".^ ~ MHz. This expected scale is about twice the ob- 
served scale Avjss- However, in |3 above we show that for 



PSR J0437-4715 the slopes of the time and frequency struc- 
ture functions are best described by the refractive model. The 
characteristic frequency scale is thus determined by the expres- 
sion (Shishov et al. i2003<) : 

Avlf ^ cliRerefOsc) (12) 

For a refractive angle of, for example, 9ref ~ 20sc at observ- 
ing frequency /** - 330 MHz we obtain for the decorrelation 
bandwidth: Av" , ^ 15 MHz. This value is consistent with the 

ref 

measured value. 

However, the value of Osc ~ 0.3 mas is much smaller than 
the value of Osc ~ 2 mas estimated from the statistical depen- 
dence of Osc on DM, given in the paper Pynzar' & Shishov 
(1997) using the dispersion measure of PSR J0437-4715, 
DM = 2.65 pc/cm^ (see Paper 1). Therefore, Component A 
explains the observational data only with difficulty. 

5.2. Model for Scattering by Component C 

We now suppose that component C of the interstellar medium 
makes the greatest contribution to the scintillation of PSR 
J0437-4715. Using a distance of R = 10 pc for compo- 
nent C, at /" - 330 MHz we obtain for the Fresnel scale 
bfr = {R/ky^^ ^ 2 X 10'" cm. Using an observer speed of 
y = 30 km/s, dominated by the Sun for this nearby material, we 
find for the characteristic spatial scale of the diffraction pattern 
bdif = Vtjss ~ 3 X lO' cm. Scintillation is again in the saturated 
regime in this case. Using this value of b^nf, we estimate the 
scattering angle as Osc ~ ^l(kbdif) ~ 1 mas. Substituting this 
value of the scattering angle in Eq.[n we obtain the estimated 
scintillation bandwidth for the diffractive model Av^.^ ^ 40 
MHz. However, as argued above, the time and frequency struc- 
ture functions favor a refractive model. For consistency with 
the measured value of decoiTelation bandwidth Av" , - 16 

ref 

MHz at /° = 330 MHz, we take 6I„/ = l.SOsc and obtain 
Av" , = 15 MHz. In this case, the value of Osc ~ 1 mas cor- 
responds much better than for the Component A model to the 
expected value of Osc ~ 2 mas for this pulsar. We therefore 
conclude that component C can explain the observational data 
better. 

5.3. Comparison of Scintillation Parameters for PSR 
J0437-4715 and PKS 0405-385 

Quasar PKS 0405-385 is located about 10° from PSR 
J0437-4715 on the sky. The interstellar scintillation of this 
quasar is certainly determined by component C of the turbulent 
interstellar plasma (Dennett-Thorpe & de Bruyn 2002 Rickett 
et al. I2OO2I 1. The scintillation of PKS 0405-385 is weak at 
/ = 4.8 GHz, with a measured scintillation index of m ^ 0.1 
(Rickett et al. 2002 1. Prokhorov et al. jl975> showed that in 
weak scintillation, the square of the modulation index is about 
equal to the value of the phase structure function at the Fresnel 
scale: 

Ds(bFr) ^0.01 for/ = 4.8 GHz. (13) 
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For a power-law spectrum with power-law index n, the value of 
Dsibfr) at frequency /o can be scaled to the given frequency / 
by the relation: 

Ds{f,bFr(f)) = {hlft^'-^'^Ds{h,bFr{m (14) 

Here we used the frequency dependence of modulation index 
(Shishov'19931. Using this relation with n - 3.6 and the above 
value for the phase structure function at / = 4.8 GHz, we ob- 
tain the estimate Dsibfr) ~ 20 at /o - 330 MHz. We can also 
scale the timescale for weak scintillation f^, = 33 min at 4.8 
GHz (Rickett et al. .2002j to fo - 330 MHz using the relations 

tcr « (fo/fcry^ (15) 
tiSS ~ (fo/.fcr)tcr, /() < fcr- 

Here /„ is the frequency that divides weak and strong scintil- 
lation regimes. It is given by the equation 

Ds(fcr,bfr(.fcr)) = (.fo/fcrf^^^'^Ds(fo,bfr(fo)) = 1. (16) 

Using this equation and the fact that Ds{f,bfrif)) = 0.01 at 
/ = 4.8 GHz, we obtain = 1 GHz, f„ = 70 min. The 
predicted time scale of tjss - 23 min at frequency fo = 330 
MHz is close to the value obtained for PSR J0437-47I5. This 
suggests that the same scattering medium causes the flux vari- 
ations of quasar PKS 0405-385 and the scintillation of PSR 
J0437-47I5. 

Another, stronger indication of this can be obtained by in- 
cluding the point corresponding to quasar scintillation param- 
eters in Fig. 12 (triangle). Here we use a time scale of 33 min 
and the value of Ds - 0.01 scaled to /o - 1 GHz. The point 
is in a good agreement with extrapolation from the pulsar data. 
A fit to all points gives the slope a - 1.8 + 0.15 (the solid line 
in Fig. 2). Using this model of the local scattering layer (with 
y = 30 km/s) we can reduce the time and frequency structure 
functions to a spatial form. Fig. |2l (top x axis) shows the in- 
ferred phase structure function, using conversion of the of time 
scale to spatial scale as b - 3 x 10^ cm s^^A/. 

6. Conclusions 

Observers have reported scintillation bandwidths for PSR 
J0437-4715 that range by a factorof about 30, when scaled to a 
single observing frequency. The observations of Paper 1 found 
two scales, at each of two observations. In this paper, we show 
that all of those observations can be represented by a single 
power-law structure function, when corrections for incomplete 
sampling in time and frequency are applied. Results of broad- 
band observations correspond to the actual scales of interstellar 
scintillation, including the characteristic scale at the observ- 
ing frequency. Results of narrow band observations correspond 
to the actual interstellar scintillations only for small frequency 
lags, Af « B. Thus, in many cases, the observational papers 
presented only part of the information carried by the structure 
function; in these cases we reconstructed the missing informa- 
tion via theoretical estimates. Scintillation bandwidth and scin- 
tillation timescale represent only two parameters from the ex- 
tensive information carried by the autocorrelation function; in- 
deed, the smaller of the two scintillation bandwidths found at 



each observing epoch in Paper 1 were obtained by restricting 
the bandwidth to a smaller range. 

Additional observations, that record modulation index at 
the central peak and the functional form of the structure func- 
tion, can test our suggestion that all scintillations result from 
a single power-law structure function. Alternatively, the wide 
range of scales might reflect an interstellar Levy flight as pro- 
posed by Boldyrev & Gwinn (2003). We respectfully request 
that observers preserve a greater portion of the rich information 
carried by the autocorrelation function of intensity in future ob- 
servations. 

Further analysis of the scintillation data for PSR 
J0437-4715 indicates that they can explained by a model for 
a layer with enhanced turbulence, at distance of about 10 pc 
from the Sun. This model appears to work better than models 
for scattering by more distant material. The flux variations of 
the extragalactic source PKS 0405-385 located nearby pulsar 
are consistent with scattering by material at the same distance 
and with the same scattering strength. 
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Appendix A: Structure Function for Limited Data 

Estimation of the structure function of intensity fluctuations 
has peculiar aspects, for observations limited in time span T 
or bandwidth B. The main problem lies in proper normaliza- 
tion of the correlation function. For observations with limited 
T and B, the correlation function calculated from observations 
is 

Ci(At,Af) = B,(At,Af)/Bo, (A.l) 
where Bq = B/(Af = 0, A/ = 0). 

However, the desired correlation function should be normalized 
by the square of the ensemble-averaged intensity: 

Kj(At, Af) = B,(At, Af )/{I)l,. (A.2) 

Here, the subscripted angular brackets {...}e„s indicate an en- 
semble average. The corresponding conversion factor is the 
square of the modulation index: 

mj = Bo/{I)l„. (A.3) 

The subscript "/" on m indicates averaging in the frequency 
domain; analogously, "f" denotes the time domain. 

For values of the time interval T, or frequency bandwidth 
B, that are large in comparison with the characteristic time and 
frequency scales of scintillation tjss and Avjss, the modula- 
tion index = 1 and Cj = Kj. For small values of (T/tiss) 
or (B/Avjss), the intensity fluctuations along these dimensions 
are intrinsically non-stationary stochastic processes, the vari- 
ance depends on the averaging interval, and the estimated scin- 
tillation index is less than one. 
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In the case where the averaging procedure coiTesponds to 
integration in the frequency domain, the variance of the inten- 
sity fluctuations is determined by the equations: 



{(My 



= df'[l{f')-{I)f\ , Where (A.4) 

</)/ = (l/B) r df'Iif). (A.5) 
Jo 



Here the angular brackets with subscript "/" denote integration 
over frequency as the averaging procedure. The variance (IA.4> 
is a random value. Forming the ensemble average, we find 



(A.6) 



1 

mif}f}e„s = J dfij df2[{l\„., - {Iifl)I(f2))e„s] 

1 

" "02 Jo Jo '^■^2[B/,.«.,(0)-B/,«„.(/i -/2)]. 

where Bj^ensifi - fi) is the intensity correlation function aver- 
aged over the ensemble. We have added the subscript "ens" to 
specifically denote the ensemble averaging. For a refractive or 
diffractive model forming the frequency structure of the difFrac- 
tive scintillation, the correlation function can be presented by 
the expression ('Shishov f2003> : 



Blensifl - fl) = {I)LsXp[-Ds(fl -fl)] 



(A.7) 



The dependence of Dsif\ - fi) on (f\ - fi) is different for re- 
fractive and diffractive models. For the diffractive model, we 
have (ShishovEOnill: 
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£>s(/i-/2) = 0.7[(/i-/2)/Av„s] 



(«-2)/2 



(A.8) 



Using this function we obtain 

(mj) = [2.8/(«(« + 2))](B/Ay„s)*'"''^'- (A.9) 
For the refractive model we have: 

Dsifi - fl) = 0.7[(/i - /2)/Av„s (A. 10) 
Using this function we obtain 

(mj) = [0.7/(«(«- 1))](B/Ay„s)<"-''. (A. 11) 

The measured value of the scintillation index is a partly 
averaged one, and in the case of small values of (T/tiss) or 
(B/Av/55) it can fluctuate, by about 100%. 

From similar considerations for the temporal version of the 
averaging procedure we obtain 

(mj) = [l/(«(«- l))](r/f„s)<"-''. (A.12) 

If we employ two-dimensional integration, in frequency 
and time, the scintillation index is determined by either Eq. lA.9l 
or lAJll for (T/tiss ) « (B/Av,ss), and by Ea. lAl2l for the re- 
verse case. 
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Fig. 1. Measurements of the scintillation bandwidth and scin- 
tillation time scale plotted against observing date. Results are 
scaled to /o = 330 MHz (left vertical scale) and /o = 1 
GHz (right scale) for purposes of comparison. The horizontal 
lines indicate the effective bandwidth used for measurements of 
Av/ss . Symbols show reference: squares: Nicastro & Johnston 
JiggS*): circles: Johnston et al. (1998); triangle: Issur (2000i; 
stars Paper 1 (4-pointed: broad scintillation, 5-pointed: nar- 



row); cross Gothoskar & Gupta J2000t . 
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Fig. 2. The time structure function of phase fluctuations for 
PSR J0437-4715 reduced to the reference frequency /o - 
1000 MHz, as compiled from the observations. Open sym- 
bols are as in Figure 1. The solid triangle indicates scintilla- 
tion of the quasar PKS 0405-385 Rickett et al. (EUU^ . The 
dashed line indicates the best fit to data for the pulsar (index 
ai = 1.46 + 0.20); the solid line corresponds to the best-fitting 
power-law to all points (index a2 - 1.8 + 0.15). Top x-axis 
corresponds to spatial scale of the inhomogeneities. 
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Fig. 3. The frequency structure function of phase fluctuations, 
reduced to the reference frequency fo = 1000 MHz, as com- 
piled from the observations, using the same symbols as in 
Figure 1. The solid lines show the diffractive model (upper) 
and the refractive model (lower). Slopes were fit to the data 
taken at/ = 327 MHz. 
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Table 1. Measurements of Scintillation Parameters for 
J0437-4715 
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" Scaled to v = 330 MHz using Ay oc / 

References: (1) Nicastro & Johnston 1995; (2) Johnston, Nicastro, & 
Koribalski 1998; (3) Issur 2000; (4) Gwinn et al. 2005 (Paper 1); (5) 
Gothoskar & Gupta 2000. 



Table 2. Expected Modulation Indices for Observations of 
PSRJ0437-4715 
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References: (1) Issur 2000; (2) Gothoskar & Gupta 2000; (3) Gwinn 
et al. 2005 (Paper 1); (4) Johnston, Nicastro, & Koribalski 1998. 



